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We calculate nuclear-polarization energy shifts for hydrogenlike ^gfPb*^''". A retarded 
transverse part as well as the Coulomb part is taken into account as the electromagnetic 

o ' 

. interaction between an electron and the nucleus. With a finite charge distribution for 

<N : 

the nuclear ground state and the random-phase approximation to describe the nuclear 

o 

' excitations, we obtain nuclear polarization energy of the lsi/2 state as -38.2 (-37.0) meV 



in 



> 

(N 



X 



in the Feynman (Coulomb) gauge. For the 2si/2, '2p\/2 and 2p3/2 states, they are -6.7 (- 
6.4), -0.2 (-0.2) and +0.0 (+0.0) meV, respectively. The seagull term in the two-photon 
exchange diagrams is shown to be quite important to obtain the gauge invariance of 



o 

' nuclear polarization energies. 



o 

^ ; PACS number(s): 31.30.Gs, 31.30.Jv, 12.20.Ds 

CJ \ I. INTRODUCTION 



High-precision Lambshift measurement on liigh-Z hydrogenlike atoms |l] arose a renewed interest in the 
quantum electrodynamic (QED) calculation of electronic atoms. Comparison of theoretical results with 
corresponding experimental data allows sensitive tests of QED in strong electromagnetic fields H,^. 
Any discrepancy between theory and experiment may either motivate an improvement of theoretical 
calculations and a refinement of experiments, or it may indicate a possible influence of non-QED effects. 
In this context, the study of nuclear polarization (NP) contributions to the total energy shift of atomic 
levels becomes important because as a background effect, it represents a natural limitation of any high- 
precision test of QED. Unfortunately, evaluation of NP is not practicable on the first principle. Any 
calculation of NP is inherently phenomenological and depends on the parameters of the nuclear model 
used to describe the intrinsic nuclear dynamics. 
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During the past years, a lot of experience has been accumulated in calculating the NP effect for muonic 
atoms 1^ . There it leads to a large correction at a keV level, mainly because of the huge overlap of the 
muon-wave function with a nucleus and because the transition energies in muonic atoms are of the order 
of magnitude of typical nuclear excitation energies. 

Much less attention has been paid to the NP effect for electronic atoms. They turn out to be reduced 
by orders of magnitude because of the small overlap of the electron- wave function with the nucleus and 
because the transition energies in electronic atoms are in general orders of magnitude smaller than typical 
nuclear excitation energies. 

The NP effect for electronic atoms was first calculated in terms of the second-order Schrodinger pertur- 
bation theory A relativistic field-theoretical treatment of nuclear polarization calculation was then 
presented by Plunien et al. |^-^ utilizing the concept of effective photon propagators with nuclear polar- 
ization insertions. The formalism allows to take into account the effect of the electron negative-energy 
intermediate states besides the usual contribution of the electron excited into higher unoccupied inter- 
mediate states. They found that in electronic atoms NP energies become small due to the cancellation 
between contributions of positive energy states and those of negative energy states. 

In the above studies, only the Coulomb interaction was considered based on the argument that the 
relative magnitude of transverse interaction is of the order of (w/c)^ and the velocity v associated with 
nuclear dynamics is mainly nonrelativistic. However, it may not be justified. In fact the importance of 
the transverse interaction has been reported for muonic atoms [^,0 . In Rcf . , the transverse nuclear 
polarization has been studied in order to explain the discrepancies between theory and experiment in 
the 2p and 3p fine-structure splitting energies of muonic ^gfPb. The contribution for the muonic lsi/2 
state amounts to 20% of that of the Coulomb interaction. 

The transverse interaction is expected to be much more important for electronic atoms than for muonic 
atoms because of its long-range nature. The transverse nuclear polarization for heavy electronic atoms 
was first studied by Yamanaka et al. 12| using the Fcynman gauge and a collective model for the 



nuclear excitations. They found that the transverse contribution is several times larger than the Coulomb 
contribution in heavy electronic atoms before the contributions of the positive and negative energy states 
cancel. However, due to nearly complete cancellation between them, the total NP energy becomes very 
small. 

The purpose of the present paper is twofold; The one is to see how much NP energies are expected 
for a practically best available model of the ^gfPb nucleus. For this purpose the Dirac-electron wave 



functions are solved in the Coulomb potential with a finite nuclear charge distribution and the random- 
phase approximation (RPA) is used to describe the nuclear excitations. The other is to see whether 
NP energies are sensitive to the choice of the gauge. For this purpose NP energies are calculated both 
in the Feynman and Coulomb gauges. We will see that NP calculation with only the ladder and cross 
diagrams shows large gauge dependence and inclusion of the seagull diagram removes most of its gauge 
dependence 

Calculations are carried out in momentum space. They involve only double integrals, which are easily 
carried out with high precision. 

II. NUCLEAR POLARIZATION CALCULATION 

The second-order contributions to the nuclear polarization are given by three Feynman diagrams in 
Fig. 1. (Here we regard the seagull graph as one of the nuclear polarization diagrams.) Two photons are 
exchanged between a bound electron and a nucleus and the nuclear vertices are understood to have no 
diagonal matrix elements for the ladder and cross diagrams, and no nuclear intermediate states for the 
seagull diagram. 

The nuclear-polarization energy shift due to the ladder and cross diagrams is given by , 
AEnp = z(47ra)2 [ d^xi ■ ■ ■ d'^Xii,{xi)-f^S'Fixi,X2h''^ix2) 



xD^^{xi,X3)Il^j^{x3, X4) Dq^{x4, X2). (1) 

Here ip is the electron wave function, 5*1^ the external-field electron propagator, D^^ the photon propa- 
gator and h|^ is the nuclear polarization tensor which contains all information of nuclear dynamics. We 
use units with h — c — 1 and — Aira. 

In terms of transition charge-current densities, the electron and nuclear parts of Eq. (|l|) are written as 

J 27r ^ E-LUe + iEi>e 



and 

doj _^-^^(t,-u)^ I Jn{x3)iI'Jn{x4)i'I JN{Xi)lI'J%{x3)ri 

I' 

where lo,, — Eii — Ei and lon = Eji — Ej are excitation energies of electron and nucleus, respectively. The 
suffixes i{I) and stand for the initial and intermediate states of the electron (nucleus), respectively. 



n^^(a;3 X4) = / —e-''^<'*^-^'^^ V ( ■^n(^3)ii'Jn{x4)i'i _ Ji{x4)ii'Ji{x3)i'i \ 
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In the momentum representation, the NP energy shifts due to the ladder and cross diagrams are 
written as 



i' I' 



and 



i' 

respectively. The substitution 



J^[-q)zt'Je{Q )^'^ J kj-l ) 1 1' -^nW I' I (5) 

i' I' 



nip 

in Eq. (|]) gives the energy correction due to the seagull diagram ||l3| , [l^ 

Je(-g)»»'Je(gO»'» PN{q-q')ii ^a 

i' ^ 

Here is the proton mass, S^'' the Kronecker delta extended to four dimensions with (5'^" — 0, and 
Pn{x)ii is the ground-state charge distribution of the nucleus. The total NP energy shift is given by the 
sum, AS^p + AE^p + AE^%. 

Substituting the corresponding multipole expansions for the Fourier transforms of the currents, these 
NP energy shifts are written in terms of the multipole form factors of the electron and nucleus defined 

by 



< i' II mx{q) \\i> = j dx ]x{qx) < i' \\ Yx{n^)pe{x) \\ i >, (8) 

< i' II txL{q) \\i> ^ Jdx ]L{qx) < i' II Yxl{^^) ■ j^{x) II I > , (9) 
for the electron, and 

< /' II Mx{q) \\I> = j dx jxiqx) < I' \\ Yxin^)pNix) \\ I >, (10) 

< /' II TxLiq) \\I> = Jdx jL{qx) < I' II YxLin^) ■ Jn{x) II / > , (11) 



for the nucleus. In the above equations, jx{qx) is a spherical Bessel function, Yxl is a vector spherical 
harmonics and A is the multipolarity of transition. With these substitutions, angular parts of q and q' as 
well as uj integrations in Eqs. (^), (^) and can be carried out analytically. Integrations with respect 
to q and q' are carried out numerically. Nuclear-polarization energy shifts are thus given by the sum of 
these double integrals over the nuclear and electron intermediate states. For the seagull contribution, 
the summation is only over the electron intermediate states. 

In the following, we shall give formula for the NP energy both in the Fcynman and Coulomb gauges. 
We restrict ourselves to = 0+ for the spin-parity of the nuclear ground state. The spin-parity of the 
nuclear intermediate state /"^ is in this case equal to the spin-parity of the transition . 

A. The Feynman gauge 

The photon propagator in the Fcynman gauge is given by 

D^^^{co,q) = -^f-. (12) 

Here the metric tensor 17^^ is defined by goo = 1 and ga = —1. With this propagator, the NP energies 
due to the ladder, cross and seagull terms are given by 

AEkp = -E (2, + i)(2/^ + i) (^)'^ X ^'^ X ^^^'^''^'^ W^iq)W^{q')5i„, (13) 
AE^P^-Y,—-^^^^^—-i-rV Hdq Hdq' I.{q,q')W^{q)W^{q')5vx. (14) 

00 pOO 



^^^p = - E #?tt(-)'^ r^^ r^^' 

V + 1) TT Jq 7o 



)>V^g(9,'z), (15) 



where V denotes the Cauchy principal value of any improper integral. The q and q' integrals may be 
improper in case of the NP calculation for the electron excited states. 
In the above expressions, /+,/- and Isg ai'c functions written as 

'^'^ " 2{Q, + q){Cj, + q%UJM + q){^N + q') ^^^^ 
Sgn{E,)[^ + (^e + + g + q')] ± ± 9{±E,)^^ 

q + q q + q lun + t^e 

T I '^ IJ7 \ qq'i^e + q + q') 

IsG{q,q ) = sgn{E,,)- ^ - ^^ - ^ ' ^^^^ 

2mp (q + q')(q + LUe){q + uJe) 

where a)e = sgn{Ei')uje, while Wj^{q), Wxiq) and yVscil) written in terms of the electron and 

nuclear form factors as 



Wf (q) = 5^[< i' II mx{q) \\ i X I' \\ Mx{q) \\ I > 

A 

A+1 



_ ^ (-1)^+1-^ < i' II txi^iq) II i >< I' II T,r.{q) \\ I >], 

L=A-1 

W| (9) = II II * X ^' II Mx{q) II / > 

A 

A+1 

+ < II iAL(g) II i X /' II II / >], 

L=A-1 

Wsg{<i,q') = E E < II IN >< ^ II ^^(9' 9') II I X i' II iAL(g') \\i>. 



A L=A-1 

The seagull term contains the Fourier transform of the nuclear ground state, 



< / II ML{q,q') \\I>= J r'^drpN{r)ii3L{qr)3L{q'r). 



B. The Coulomb gauge 



The photon propagator in the Coulomb gauge is given by 



and {5ij — in Dfj projects out both transverse parts of electronic and nuclear currents. 
Making use of the relation 

we obtain the NP energies for the ladder, cross and seagull terms: 

+ I+{q,q')W?{q)W?{q') + I^^ {q')WE {q)W? {q') } ^/'A, 

+ I-{q,q')W?{q)W?{q')+I^'^{q')WE{q)W?{q')}Srx, 
AEl,% = -J2j^drV Hdq rdq'IsG{q,q')W'iGM). 



In the above expressions, lY^{q') and ib^ {q') come from the interference between the longitudinal 
and transverse contributions and given by 



tLT, _ ,sgn{E,,)q'iu, + lon) ± e{±E,,)2q'^ 



{q' + U)e)iq' + t^Jv)(We + t^jv) 

In the Coulomb gauge, W^{q), W^{q) and W^{q) are given by 

= E < II "^^(«) II * >< ^' II ^'^^('?) II ^ >' 



(27) 



(28) 



A+l 



E <^'IUAL(<7)|h></'|lT,i(g)|l/> 



L = A-1 



WFg(9,9') = E 



E (< II UxL{q) II * X / II MLiq,q') \\ I >< || ^iAL(9') || * >) 



A LL=A±1 



(29) 



+ <i'\\ txxiq) II I >< I II Mx{q,q') \\ I >< i' \\ txxiq') \\ i > 



with 



< i' II uxx-i{q) II t > =<i' II tx-i{q) || i > -i 



A UJe 

2X+1 q 



< i' II mx{q) II i >, 



< II "AA+i(g) II j > - < II t\+i{q) II « > < II ™^('^) II * > 



(30) 

(31) 
(32) 



C. Electron wave functions 



The radial Dirac equations for the electron are written as 



ar r ' 



(33) 



d K 



Fe,. = [m,-E + Vir)]GE.Ar), 



(34) 



dr r ^ 

where the potential V{r) is obtained from the the ground-state charge distribution of ^''^Pb, which is 
assumed to be a two parameter Fermi distribution 



PN{r) 



II 



1 + exp[(r - Ro)/c 



(35) 



with Rq = 6.6477 fm and a = 0.5234 fm |10|. These equations are solved numerically by using the 
fourth-order Runge-Kutta method. For both positive and negative energy continuum states, the radial 
functions are normalized by the conditions 



1/2 



r — >oo \ Tip 



sin{pr + 5) 



(36) 




1/2 



cos(pr + S), 



(37) 



with p — — ml . The bound-state wave functions are normahzed as [/g°°(G^^ + F^^^)dr 1]. 

Transition form factors for the electron, (H) and (||), are calculated using the formula given in Ref. |]l0|] . 
They are stored in the computer with six different step sizes of Aq depending on the electron energy 
Eii . In Fig. 2 we show the E\ charge form factors < £'i',pi/2 || mi{q) \\ lsi/2 > of the electron with 
three different energies, Eii—2, 6 and 10 MeV. One finds that they have sharp peaks at g = Eii and 
decrease rapidly as q increases. Numerical integrations in Eqs. (p^)-(p^) and (p4[)-([2^) are performed 
by Simpson's one-third rule. 

Most of the NP correction comes from the continuum states with energies greater than rrie and energies 
less than —me. Summation over the electron states i' in Eqs. (p^)-(p^ and Eqs. ([2^)-(|2^) implies an 
integration with respect to Ei' . The Eii integration is carried out by using the Gauss-Legendre quadrature 
over the intervals -250MeV< E^i < —me and m^ < Eii <250MeV. Accuracy of the numerical results 
is checked by comparing them with the results of Simpson's rule. We also included the electron bound 
states for n'<7 in the calculation. 



The random-phase approximation is used to describe the nuclear excitations. The calculation we 
employed is the same as those performed earlier in Ref. i.e., the same single-particle basis, the same 
particle-hole configuration of approximately a full 3huj space and the same Migdal force p5| ] parameters 
to describe nuclear two-body interaction. Nuclear transition form factors , ( [l0| ) and (|Tl]), are calculated 



The calculated charge and magnetic current densities are examined by comparing with experimental 
B{E\) and B{M\) using the following relations; 



D. RPA calculation of nuclear charge and current densities 



by using the formula given in Ref. . 




(38) 



and 




(39) 
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In Table I, we compared the energy-weighted sum of B{EX) over the RPA states with the classical 
energy- weighted sum-rule value (EWSR) of Ref. jl6|. The results for the £^0 and El transitions exceed 
the EWSR by 20% and 10%, respectively while the results for the E2 and E3 transitions agree well with 
the EWSR. For the Ei and E5 transitions, our results exhaust only 50% of the EWSR. This may be 
due to the insufficient configuration space for the EA and E5 calculations. 

On the other hand, there is no experimental constraint imposed on the nuclear electric current. How- 
ever, the electromagnetic current should satisfy the continuity equation required by charge conservation; 

^Jn + VJn = 0. (40) 
ot 

Using the nuclear Hamiltonian H^^ we rewrite it as 

i[HN,PN]+yJN ^0- (41) 

Taking the matrix element of Eq. ( pl| ) between the initial and the final nuclear states, we obtain charge 
conservation condition; 

Here, lon = Eji — Ej is the energy of nuclear excitation. 

For any kind of model calculation involving the nuclear current, it is necessary for the model to 
satisfy the charge conservation condition in order to observe the gauge invariance p3[ |. Unfortunately, 
the charge-current densities constructed from the present RPA calculation do not satisfy the charge 
conservation of Eq. ( ^2| ) . The violation of the charge conservation comes from the inconsistency of using 
empirical single-particle energies together with the impulse charge-current operators, as is discussed in 
Refs. [ P^IJTtI . It is desirable if one could construct a microscopic self-consistent model together with the 
nuclear current satisfying the charge conservation which is realistic enough to reproduce the observed 
spectra and B{EX) values. However, the refinement of calculation will be left for a future work and at 
present we are satisfied with the fact that the calculated NP energies show only a small gauge violation 
even though the empirical single-particle energies are used in the RPA calculation. 

III. RESULTS AND DISCUSSION 

NP energy shifts are obtained by computing an energy shift for each of the RPA excitations and 
summing the results. Our calculation gives 38, 129, 160, 222, 202, 218 and 70 nuclear states for the 0+, 
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1^, 2+, 3^, 4+, 5^ and 1+ excitations, respectively. Figure 3 shows the NP energy spectra of the lsi/2 
state for the respective nuclear spin-parities. These spectra are calculated in the Coulomb gauge with 
the Coulomb and transverse parts of the electromagnetic interaction. They are very similar to the RPA 
spectra of B{EX) and B{Ml). 

Tabic II summarizes the NP energies of the lsi/2 state for ^g|Pb^^+. The 1st column denotes nuclear 
spin-parities. The entries in the 2nd column indicate the contributions to the NP energy from the ladder, 
cross and seagull terms as well as those of the positive and negative energy intermediate-states of the 
electron. The 3rd column shows the NP energies in the Feynman gauge, while the 4th column shows 
the NP energies in the Coulomb gauge. The transverse contributions are included in both columns. The 
5th column shows the Coulomb NP energies without the transverse contribution (hereafter referred to as 
CNP). The 6th column shows the results of the previous NP calculation in the Feynman gauge assuming 
a collective model for the nuclear excitations [|l^ . The 7th column is the CNP of the same calculation. 
Finally, the last column shows the CNP calculated by another group [^. 

In Table II, it should be noted that the positive and negative energy continuum states of the electron 
always contribute to a NP energy with opposite sign, the fact that is already observed in Rcfs. ||^,p^. 
It is also understood from the fact that the contribution of the negative energy electron describes the 
blocking of the response of vacuum due to the occupied state of the atomic electron. In electronic atoms 
the virtual pair creation requires energy of only 2me, which is smaller than the typical excitation energy 
of a nucleus and is correspondingly important. 

The 5th column (the present CNP) may be compared with the 7th and 8th columns of the previous 
calculations. All three calculations show very good agreement with one another. The agreement indicates 
that provided the B{EX) values are similarly chosen, CNP of electronic atoms is not very sensitive 
to the detail of transition charge densities. The monopole NP energy shows some difference between 
the calculations. It is -4.0 meV (5th column), -7.2 meV (7th column) and -3.3 meV (8th column), 
respectively. The difference is mainly caused by the Dirac-electron wave functions used; those for the 
finite charge distribution are used in the 5th and 8th columns, while those for a point charge are used 
in the 7th column. The difference is conspicuous only for the monopole NP energy because the nuclear 
monopole transition-potential exists only inside the nucleus where the electron- wave functions generated 
from the point charge and the finite charge distribution differ appreciably. 

The dipole NP energy also shows some difference, i.e., it is -20.3 meV (the 4th column), -19.5 meV 
(the 6th column) and -17.6 meV (the 7th column), respectively. The difference between the 4th and the 
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7th columns arises from the fact that the energy-weighted sum of B{E1) over the RPA states exceeds 
the classical EWSR value by 10% (see Table I). 

The NP energies including the transverse effect (the 3rd column) can be compared with the results of 
Ref. ijl^ (the 6th column), both of which were calculated in the Fcynman gauge. Agreement between the 
two calculations is good except for the nuclear monopole excitation. Here again, we may conclude that 
NP energies of heavy electronic atoms are not very sensitive to the detail of transition current densities. 

An important feature of the present calculation, which is in fact crucial for the numerical estimate of 
NP energies, is that there exists a large violation of gauge invariance in the NP energy-shifts as far as 
only the ladder and cross diagrams of Figs. 1(a) and 1(b) are taken into account. By using the minimal 
prescription, the nonrclativistic electromagnetic interaction involves the square of vector potential called 
a seagull term coming from the kinetic energy, and this term is necessary for the gauge invariance in a 
nonrclativistic system |]l3| , |l4| . It is interesting to investigate whether the inclusion of the seagull term 
restores the gauge invariance of the present RPA calculation. (A proof of the gauge invariance is given 
in the Appendix.) 

As is seen from Table II, this is nicely confirmed numerically. With the ladder and cross diagrams, the 
NP energy for the lsi/2 was +1.5 meV in the Feynman gauge, while -32.7 meV in the Coulomb gauge. 
The gauge dependence was 34.2 meV. However, after the inclusion of the seagull term, it is -38.2 meV 
in the Fcynman gauge and -37.0 meV in the Coulomb gauge. The gauge dependence is reduced to 1.2 
meV. This small gauge dependence shows that the seagull term is quite important in restoring the gauge 
invariance of the NP calculation. The fact also implies that the use of empirical single-particle energies 
in the RPA calculation does not introduce a serious violation of gauge invariance into the NP energies. 

In Table II, we recognize that the nuclear dipole states give predominant NP contributions. The 
transverse contribution in particular is very large if the contribution from the positive and negative 
energy states are separately considered. For example, in the Coulomb gauge, the Coulomb plus transverse 
contribution from the ladder and cross diagrams is -120.5 meV, while the corresponding contribution 
without the transverse term is -37.0 meV. The reason for this fact is understood as follows. For simplicity, 
we disregard the retardation of the interaction. Then Eq. (|2^ ) gives the Coulomb-Breit potentials for 
electron: 




(43) 



(44) 
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In Fig. 4 the radial parts of the potentials due to the 14.6 MeV state are plotted after multipole 
expansions of Eqs. (|4^ ) and (Q). The Coulomb potential is larger than both of the transverse potentials 
in the region r < 60 fm, while the transverse potential with A = 1,L = becomes larger than the Coulomb 
potential in the region r > 60 fm. This long-range nature of the transverse potential explains the fact that 
the El transverse contribution is dominant for electronic atoms because the Bohr radius of the lsi/2 state 
is much larger than 60 fm. Since the radial dependence of the Coulomb potential (transverse potential) 
with multipolarity A is l/r^^^ i^/f^) outside the nuclear region, the magnitude of the potentials in the 
region of the Bohr radius decreases as the multipolarity increases. Therefore, for multipolcs other than 
the dipole, i.e., A'^=2+, 3^, 4+ and 5^, the NP energy is mainly determined by the overlap between the 
electron transition density and nuclear multipole potential in the region of the nuclear surface, hence 
the transverse contribution is small compared with the Coulomb contribution. The fact also explains 
that the Coulomb contribution is predominant for heavy muonic atoms in any multipolarities because 
the overlap near the nucleus is always significant due to the small muon Bohr-radius . 

The dominance of the El contribution in the NP energy can be seen more clearly in the spectral 
density of the NP contribution from a particular nuclear excitation as a function of electron energy. 
In Fig. 5, these spectral functions in the Coulomb gauge are shown for three different nuclear states; 
(a) 0+ (13.3 MeV), (b) 1" (14.6 MeV) and (c) 2+ (10.2 MeV). In each panel, the solid line shows the 
spectral function including the transverse contribution, while the dotted line shows the result without 
the transverse effect. The NP energy due to each of the nuclear states in Fig. 3 is given by the integral 
of the corresponding spectral function over the electron energy. One can see that the low energy region 
of the spectral function for the 1^ state (Fig. 5(b)) is different from the other two. The El spectrum 
shows a peak at threshold {Eii — me). It should be reminded that only the < || Jio(g) |j / > and 
< i' II ^10(9) II * > a-re nonvanishing at g = 0, so that the large overlap between these transverse form 
factors is guaranteed in the low momentum region. 

For the excited L-shell electrons, we can repeat the discussion that the transverse El multipole plays a 
crucial role in the NP effects of hydrogenlike atoms and essentially determines the magnitude of the NP 
effects. The total NP shifts for the lsi/2, 2si/2, and 2^3/2 states are summarized in Table III and 

are compared with the Coulomb NP energies. The relative importance of the transverse contribution is 
about 10-20% effect for the lsi/2 and 2si/2 states, and the level shifts for the 2pi/2 and 2^3/2 states are 
negligible. 

Since the nuclear dipole states have predominant NP contributions, we must note here effects of the 
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spurious ccntcr-of-mass motion of the nucleus on the NP energies. The present Migdal force brings down 
the lowest 1^ state to the imaginary eigenvalue of 1.32i MeV. Since this 1^ state carries most of the 
spurious center-of-mass motion, we excluded this nuclear state from the NP calculation. The 0.7% of 
the spurious center-of-mass motion remains in the rest of the 1^ states, whose effects on the NP energies 
are negligible. Thus our results for NP energy due to the cross and ladder diagrams contain intrinsic 
excitations only. On the other hand, the seagull contributions calculated by using Eq. (|2^) contain both 
of the effects of intrinsic excitations and center-of-mass motion. The seagull contribution coming from 
the center-of-mass motion must be eliminated for the dipole mode. This was achieved by using the 
effective dipole charges Cp — N/A, e„ = —Z/A instead of using the true charges Cp — 1, e„ = 0. 

IV. SUMMARY 

We have calculated the NP energy shifts for the hydrogenlike ^gfPb^^"'' taking into account the effects 
of the electron in the negative energy continuum besides the usual contributions of the electron excited 
into higher unoccupied orbitals. The evaluation of the NP energies contains the seagull graph as well as 
the ladder and cross diagrams. The Dirac-electron wave functions were solved in the Coulomb potential 
with a finite charge distribution for the nuclear ground state and the RPA wave functions were employed 
for the nuclear excited states. 

The results presented in the previous section can be summarized as follows. l)In the Feynman 
(Coulomb) gauge, we obtained the NP energies of -38.2(-37.0), -6.7(-6.4), -0.2(-0.2) and -I-O.O(-I-O.O) 
meV for the lsi/2, 2si/2, 2px/2 ^-nd 2p^/2 states, respectively. 2) The nuclear dipole states have predom- 
inant contributions to the NP energy of heavy electronic atoms. In particular the exchange of transverse 
photon plays a crucial role. The net El contribution, however, becomes quite small due to the cancel- 
lation between the contributions of the electron in the positive and negative energy intermediate states. 
3) Due to the cancellation, the transverse contribution to the total NP energy shift is about 10% of 
the Coulomb contribution for the lsi/2 state. 4) The NP shifts of electronic atoms have serious gauge 
dependence if one calculates them with only the ladder and cross diagrams of the two-photon exchange 
processes. Indeed the NP energies for the lsi/2 state due to these diagrams are -1-1.5 meV and -32.7 
meV in the Feynman and Coulomb gauges, respectively. The inclusion of the seagull graph gives the NP 
energies of -38.2 meV and -37.0 meV. The seagull graph is quite important in restoring the gauge invari- 
ance of the NP calculation. 5) The NP energy for electronic atoms is not sensitive to the details of the 
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transition charge and current densities because it is primarily given by the overlap between an electron 
transition density and a nuclear multipole potential in the region outside the nucleus. 6) The present 
RPA calculation using the empirical single-particle energies together with the impulse charge-current 
does not seem to introduce serious gauge dependence into the NP calculation. 
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APPENDIX A: GAUGE INVARIANCE 



To prove the gauge invariance of the NP energy, we write the sum of the ladder and cross contributions 
as follows, 

AEh. + AE^ ^(4"")' /d- / f dq' 



-'NP ^ ^^NP 



2ti J (27r)3 J (27r)3 

N ^ 



nr (c^, q, q')D,,duo, g)i?c-(c^, g')n«i(w, 9, (Al) 



where 11^'^ (w, q, q') and 11^ (w, q, q') are the electronic and nuclear polarization tensors defined by 

nr(c., g, - E ( - m^^Mziti) , (A2) 

n«(^, q, ci)^y( JN{Q)n'Ji{-q')ri 4(-gO//'4(g)/'/ \ (A3) 

Photon propagators in the Feynman and Coulomb gauges are related to each other by 
D^Al^o.) = Dl,i,M - ^ q,q,-^iq,9,. + q,9,.) ^ 

If both g^n^"^ = and 95 are satisfied, it is easy to see that the Feynman and Coulomb gauges 



give the same result for the NP contributions given by Eq. (Al). 

Multiplying both sides of the electron polarization tensor by q^, and using the continuity equation of 
the charge conservation, one obtains 

g^nr(c^,g,g') = E(< ^\Pe{-q)\i' >< ^'\fe{q% >-< AfeiQ'W >< ^'\M-q)\i >) 

i' 

= <»i[pe(-9),i:(g')]N>- (A5) 

In deriving the second equality, we have assumed the completeness of the intermediate states of the 
electron. For the electromagnetic charge and current operators 

Pei-q) = (11) e'"-^. je(g') = o) e-^'^'-, (A6) 



used with the Dirac-electron wave functions, the commutation relation in Eq. (A5) vanishes. Hence the 



gauge invariance qfjIl^'^{LU, q, q') = for the electronic polarization tensor follows. 



For the nuclear polarization tensor, we can obtain the similar form to Eq. ( A5) by assuming the charge 
conservation as well as the completeness relation. In the present calculation, the impulse charge and 
current operators 

i i ^ i 
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are employed with the nonrelativistic RPA calculation. The spin-current operator in the second term of 
J Nil) commutes with p^. Hence the spin current introduces no gauge violation into the NP calculation 



of Eq. (Al). The convection-current, on the other hand, does not commute with pjv leading to a violation 



of gauge invariance: 



< mN{'q),JNiQ'm —PN{q-q')lI. (A8) 

mp 



Therefore the NP contribution given by Eq. (Al) is not gauge invariant with the impulse charge-current 



operators. When the seagull tensor of Eq. (g) is added to the nuclear polarization tensor of Eq. (A3), 
this term is just canceled. Hence the gauge invariance of NP calculation is restored by the seagull term 
together with the the ladder and cross terms. 
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TABLE I. The energy-weighted sums of B{EX) over the RPA states. The classical EWSR values 
are also shown for comparison. The values are given in units of e^b'^-MeV 



EA Present calculation Classical EWSR " 



'W' L97 Tel 

EV 8.15 7.38 

E2 22.2 20.5 

E3 24.4 23.4 

E4: 14.2 23.6 

E5 11.3 23.1 

Ml'^ 294 



" The radial moments < r >p in the classical EWSR are calculated by the Fermi charge distribution 

(!)• 

The £'0 operator is defined as 0{EO)^J2pr'^ / 
" The El operator is defined as 0{E1)=J2^ -l/2T3rYi^ . 
The value given in units of /in- MeV. 
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TABLE II. Nuclear polarization correction (meV) of the lsi/2 state in ^g|Pb^^+. Energy shifts AE^ , 

AE^ and AE'^'^ are contributions of the ladder, cross and seagull terms, respectively, while Ai?+ {AE ) 
denotes contribution from the positive (negative) energy intermediate states of the electron. 

present'' present^ presenf^ ReT^pP RefT^pp Ref. 
A'^ Contribution Feynman(NP) Coulomb(NP) CNP NP CNP CNP 



AE^+ 




-5.7 


-6.5 


-7.0 


AE^- 




+0.4 


+0.2 




AE''+ 




-1.2 


-0.2 




AE""- 




+2.7 


+2.7 


+3.0 


AE^+AE^ 




-3.8 


-3.9 


-4.0 






+0.7 


0.0 




AE^^^- 




-0.9 


0.0 




AE^+AE^- 


\-AE^G 


-3.9 


-3.9 




AE^+ 




-119.1 


-91.1 


-37.0 


AE^- 




+74.0 


+37.6 




AE^+ 




-49.8 


-29.4 




AE""- 




+ 110.1 


+64.2 


+ 16.7 


A 1 — 1 T A 1 — 1 \^ 

AE^+AE^ 




+15.2 


-18.7 


-20.3 






+ 144.2 


+87.8 




AE^^- 




-186.5 


-95.1 




AE^+AE^^ 


^AE'^G 


-27.1 


-26.0 




AE^+ 




-13.0 


-15.3 


-14.4 


AE^- 




+2.0 


+0.4 




AE^+ 




-3.2 


-0.5 




AE""- 




+8.1 


+9.1 


+8.6 


AE^+AE^ 




-6.1 


-6.3 


-5.8 


AESG+ 




+4.4 


+2.7 




AE^G- 




-4.1 


-2.1 




AE^+AE^^ 


^AE^G 


-5.7 


-5.7 




AE^+ 




-5.0 


-6.5 


-6.3 


AE^- 




+1.0 


+0.1 




AE^+ 




-1.6 


-0.1 




AE""- 




+3.2 


+4.1 


+4.0 


AE^ + AE^ 




-2.4 


-2.4 


-2.3 






+1.0 


+0.6 




aeS(^- 




-0.8 


-0.4 




AE^+AE^^ 




-2.2 


-2.2 




AE^+ 




-1.1 


-1.5 


-1.4 


AE^- 




+0.2 


0.0 




AE''+ 




-0.4 


0.0 




AE""- 




+0.7 


+0.8 


+0.8 


AE^+AE^ 




-0.7 


-0.7 


-0.6 






+0.5 


+0.4 








-0.4 


-0.3 




AE^+AE^~ 


^AE^G 


-0.6 


-0.6 





-6.6 



-7.2 



-3.3 



-16.3 



-19.5 



-17.6 



-7.0 



-6.3 



-5.8 



-2.9 



-2.6 



-2.6 



19 



TABLE II. (Continued). 

present" present^ present^ Ref. H'* RefTlpp Ref. 
A'" Contribution Feynman(NP) Coulomb(NP) CNP NP CNP CNP 







-0.4 


-0.5 


AE^- 




+0.1 


0.0 






-0.1 


0.0 






+0.2 


+0.2 






-0.2 


-0.2 






+0.3 


+0.1 






-0.2 


-0.1 


AE^+AE^~ 


hAS^G 


-0.2 


-0.2 


AE^+ 




-0.4 


-0.4 


AE^- 




+0.1 


+0.1 


AS^+ 




-0.2 


-0.2 


AS^- 




+0.1 


+0.1 


AE^+AE^ 




-0.4 


-0.4 






+3.8 


+3.8 


AS'^G- 




-1.7 


-1.7 


AE^+AE^~ 


^AE^G 


+1.7 


+ 1.7 



Total 



AE^+ 


-144.7 


-121.8 


-66.5 


AE^- 


+77.8 


+38.3 




A^;^+ 


-56.7 


-30.4 




AE""- 


+ 125.1 


+81.2 


+33.3 


AE^+AE^ 


+1.5 


-32.7 


-33.2 


AESG+ 


+ 154.9 


+95.4 




Ai;^G- 


-194.6 


-99.7 




AE^ + AE^ + AE^^ 


-38.2 


-37.0 





° The NP energies in the Feynman gauge. 
^ The NP energies in the Coulomb gauge. 
The unretarded NP energies in the Coulomb gauge. 

The NP energies evaluated in the Feynman gauge. Electron wave functions were solved by assuming 
the point charge for the nuclear ground state. Nuclear transition-charge densities were determined 
by a collective model. They were normalized to the observed B{EX) for the low- lying nuclear states 
and EWSR values for the high-lying giant resonances. Nuclear current densities J^w-iij") were 
obtained by solving the equation of the charge conservation ( ^ ) assuming Jnxx+i{''') = 0. 

Same as d except for the unretarded NP energies. 

■f The unretarded NP energies. Same as e except for the electron wave functions solved by assuming a 
finite charge distribution for the nuclear ground state. 
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TABLE III. Total nuclear polarization (meV) of the lsi/2, 2si/2, 2pi/2, and 2^3/2 states of ^^|Pb^^+ 
both in the Feynman and Coulomb gauges. The abbreviation CNP denotes the unretarded results. 



States 


Feynman(NP) 


Coulomb(NP) 


CNP 


lSl/2 


-38.2 


-37.0 


-33.2 


2Si/2 


-6.7 


-6.4 


-5.7 


2Pl/2 


-0.2 


-0.2 


-0.6 


2P3/2 


+0.0 


+0.0 


-0.0 
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Figure captions 



FIG. 1. Diagrams contributing to nuclear polarization in lowest order; (a) ladder diagram, (b) cross 
diagram, (c) seagull diagram. 

FIG. 2. Electronic Coulomb form- factors < Eii,pi/2 || m\{q) || lsi/2 > for the Eii=2, 6 and 10 MeV 
states (solid line). The dotted line is the nuclear Coulomb form factor for the 14.6 MeV 1^ state. 

FIG. 3. Nuclear polarization spectra as functions of nuclear excitation energy. The Coulomb gauge 
was assumed. 

FIG. 4. The transition potentials due to the nuclear 1^, 14.6 McV state. The Coulomb-Brcit propaga- 
tor was assumed. The solid, dotted and dash-dotted lines denote the potentials corresponding to pi(r), 
Jio(r) and J^(r), respectively. The transverse potential Jio(i') is seen to be larger than the Coulomb 
potential in the region r > 60 fm. 

FIG. 5. Spectral densities of nuclear polarization for (a) 0+, 13.3 McV, (b) 1^, 14.6 McV and (c) 
2+, 10.2 MeV states as functions of electron energy Eii . The solid line denotes calculation with both 
the Coulomb and transverse parts of the electromagnetic interaction, while the dotted line denotes 
calculation with only the Coulomb part of the interaction. Electron intermediate states, (a) \Eii,si/2 >, 
(b) |£'i',Pi/2 > and (c) \Ei/,d3f2 > are assumed in the respective panels. 
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